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PREFACE 


This  Train-Collision  Model  Users  Manual  presents  documentation  of  a 
simplified  mathematical  model  and  computer  program  of  the  longitudinal 
dynamics  of  rail  cars  involved  in  collisions  and  an  assessment  of  the 
range  of  subsequent  injuries  sustained  by  passengers.  The  Train-Collision 
Model,  was  originally  developed  in  the  support  of  a program  of  assessment 
of  the  crashworthiness  of  existing  urban  rail  vehicles  by  the  Calspan 
Corporation  for  the  U.S.  Department  of  Transportation,  Transportation 
Systems  Center  under  Contract  No.  DOT-TSC-681.  Further  improvements  in 
the  model  and  documentation  were  made  as  reported  herein  under  that 


contract. 
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1 . INTRODUCTION 


The  Train  Collision  Model  is  a simplified  dynamic  simulation  of 
rail-car  collisions  with  automatic  estimation  of  the  degree  of  severity 
sustained  by  passengers  of  each  car  in  a train  during  second  collisions  with 
various  types  of  interior  surfaces  at  different  initial  distances.  The  rail 
cars  are  modeled  as  rigid  bodies  separated  by  nonlinear,  energy -absorbing 
force  deflection  characteristics  with  the  first  car  of  a train  impacting  a 
rigid  barrier.  Two  modes  of  operation  are  possible,  the  first  assuming  that 
an  individual  force-deflection  characteristic  acts  between  each  car.  For 
this  mode,  the  deflection  of  each  car  (except  the  first)  is  assumed  to  be 
one-half  the  relative  displacement  between  the  two  adjacent  car  centers  of 
gravity.  The  second  mode  of  operation  allows  the  individual  deflections  of 
each  end  of  each  car  to  be  obtained  through  the  use  of  an  inertial  mass 
between  the  springs  representing  the  car  ends.  This  mode  is  useful  in 
approximating  the  different  deformations  of  cars  that  over-/underride. 

Passenger  injury-severity  distributions  are  approximated  by 
determining  the  relative  velocity  between  a mass  representing  a passenger 
and  the  car  at  times  at  which  the  passenger  has  moved  a number  of  input 
distances  through  the  car.  The  relative  velocities  are  subsequently  used 
to  calculate  a severity  index  assuming  that  a passenger  in  that  situation 
strikes  a car  fixed  object  that  nulls  the  relative  velocity  in  a number  of 
input  distances  at  a constant  force  level. 

It  must  be  recognized  that  the  mathematical  model  and  computer 
simulation  described  in  this  report  is  a very  simplified  approximation  of 
an  extremely  complex  phenomenon.  In  this  respect,  care  should  be 
exercised  in  interpretation  of  model  results  as  they  apply  to  the  real  world. 
While  it  is  expected  that  generalized  trends  in  overall  rail-car  crash 
situations  may  be  properly  obtained,  results  in  an  absolute  sense  for 
specific  crash  situations  should  not  be  indescriminantly  interpreted. 
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2.  MATHEMATICAL -MOD EL  DESCRIPTION 


2 . 1 RAIL-CAR  MODEL 


A simplified  schematic  diagram  of  the  mode  1 rail  car  model 
is  shown  in  Figure  2-1.  The  model  consists  of  a barrier  with  optional  movement 
at  a constant  velocity,  and  a number  of  rail  cars  (maximum  of  20)  idealized 
as  rigid  bodies  separated  by  nonlinear,  energy-absorbing  springs.  A 
coordinate  system  is  fixed  at  the  center  of  gravity  of  each  car  and  passenger 
at  the  instant  of  impact  and  by  inspection,  the  acceleration  of  a car  1 
relative  to  that  coordinate  system  is  given  by: 

\ = k ('A+i  ' F0  ' a=° g’ 

where  ac is  an  arbitrary  input  deceleration  (because  of,  for 
example,  train  braking). 

The  forces  acting  between  cars  arise  from  deformation  of  the  car 
ends.  Tabular  force-deflection  characteristics  are  input  to  the  model,  and 
forces  are  generated  based  on  the  relative  positions  and  velocities  of 
adjoining  cars.  In  this  mode  of  operation>it  is  assumed  that  the  deflection 
of  the  nonlinear  springs  is  shared  equally  between  cars  although  two  force- 
deflection  tables  are  input  so  that  different  characteristics  may  act  between 
different  cars.  Thus,  for  the  purpose  of  determining  the  force  acting 
between  cars  i and  i + 1,  the  deflection  and  deflection  velocity  are 


2-1 


*w 

V 


Xc  (1) 

vc  (1) 


X_  (2) 

v:«» 


X.  (3) 

v> 


BARRIER  FIRST  SECOND  THIRD 

(WITH  OPTIONAL  RAILCAR  RAILCAR  RAILCAR 

MOVEMENT) 


*p  = PASSENGER  (OCCUPANT) 


Figure  2-1  Simplified  Schematic  Diagram  of  Typical  Mode  1 Train  Configuration 
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For  the  case  of  the  first  car  of  the  train, 


h - XC1  - ^ • 


Knowing  the  deflections  and  relative  velocities  of  the  cars,  the 
forces  are  determined  as  described  in  Section  2.  2, 


A simplified  schematic  of  the  mode  2 rail  car  model  is  shown  in 
Figure  2-2.  This  mode  of  operation  allows  independent  specification  of  the 
force-deflection  properties  of  each  adjacent  car.  An  inertial  mass  is  used 
to  separate  the  cars  so  that  the  deflection  of  each  car  is  calculated 
independently.  This  mode  is  useful  in  determining  the  deflections  of 
adjacent  car  ends  when  override  occurs.  In  this  mode,  every  other  mass 
is  a representation  of  a rail  car,  thus  a maximum  of  10  cars  can  be 
simulated.  The  deflection  and  deflection  velocity  of  each  of  the  nonlinear 
springs  are  given  by: 

i+1  ci+l  ci 


i+1 


= X, 


- X, 


'i+1 


and,  for  the  first  spring, 
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Figure  2-2  Simplified  Schematic  Diagram  of  Typical  Mode  2 Train  Configuration 


2.2  CALCULATION  OF  FORCES 


The  calculation  of  forces  acting  between  rail  cars  is  based  on 
the  known  deflection  and  time  rate  of  change  of  deflection.  A general  table 
of  force-deflection  characteristics  is  used  to  interpolate  the  force  level 
when  the  deflection  is  increasing.  Unloading  and  reloading  after  unloading 
has  occurred  takes  place  along  a straight  line  with  a slope  of  URATE  as 
illustrated  in  Figure  2-3.  If  reloading  back  to  Fm,  6m  occurs,  calculation 
of  the  force  along  the  original  curve  is  resumed. 


Figure  2-3  Force- Deflection  Characteristic 


A conceptual  flow  diagram  of  the  force  calculation  logic  is  shown 
in  Figure  2-4. 
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Figure  2-4-  Force-Calculation  Logic 


2-6 


2.3  RAIL-CAR  INTEGRATION 


Integration  of  the  differential  equations  of  motion  of  the  rail 
cars  is  by  the  Improved  Euler  Method  as  illustrated  by  Figure  2-5.  The 
procedure  is  as  follows: 


At  the  beginning  of  a time  increment,  t = t , 

determine  the  acceleration  of  each  car,a^^  , 

m 

by  calculation  of  the  forces. 

Assuming  a constant  acceleration,  a^)  , over 

the  time  increment,  At  , approximate  the  velocity 
and  position  of  the  car  at  the  end  of  the  interval 
( t = tm  + At  = tm+1  ) by: 

= V + a^  At, 
m+1  m m 


= X + V At  + i a(1)  At2. 
m+1  mm  2 m 

Determine  the  acceleration  at  the  end  of  the  interval 


based  on  the  estimated  position  and  velocity,  a 
by  calculation  of  forces. 


(0) 

m+1  ’ 


Assume  that  the  acceleration  over  the  interval  is  the 
average  of  the  above  calculated  accelerations, 

and  compute  the  revised 


a = 1/2  (a(1)  + a(o)) 
m m m+1 ' * 


position  and  velocity  at  the  end  of  the  interval: 

= V + a At 
m+1  m m 

= V + j (a(1)  + a(°})  At, 
m 2 m m+1 

1 2 

= X + V At  + T a At 


m+1 


m 


= X_  + V_At  + \ (a^1}  + 
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VELOCITY 
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Figure  2-5  Improved  Euler  Method  Integration 
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2.4  PASSENGER-INJURY  ASSESSMENT 


The  severity  of  a given  crash  situation  depends  in  a general 


sense  on  the  deceleration  level  acting  on  the  occupant.  The  deceleration 
level  is  a function  of  the  type  of  surface  struck  and  the  relative  velocity  at 
impact.  Rail  cars  are,  in  general,  a relatively  uncontrolled  environment 
for  the  occupants  after  an  impact  and  thus  a deterministic  approach  to 
modeling  occupant  dynamics  is  not  at  this  time  warranted. 


severity  in  this  model  was  to  determine  a range  of  severity  based  on  the 
distance  the  passenger  travels  before  impact  and  the  crush  distance  of  the 
struck  surface.  One  passenger  is  simulated  in  each  car  and  the  position 
and  velocity  of  the  passenger  are  determined  assuming  that  an  input  constant 
deceleration  acts  (simulating  for  instance,  sliding  friction).  When  the 
position  of  the  passenger  relative  to  the  car  reaches  any  one  of  a number  of 
input  spacings  (the  distance  the  passenger  travels  before  striking  an  object), 
the  relative  velocity  is  used  to  determine  an  injury-severity  index.  It  is 
further  assumed  that  the  passenger  decelerates  to  a relative  velocity  of  zero 
in  a number  of  input  stopping  distances  (crush  distances)  at  a constant 
deceleration  level.  Note  that  the  effects  of  ridedown  (car  deceleration 
during  passenger  impact)  are  neglected  in  this  simplified  evaluation. 


The  approach  taken  to  obtain  a general  assessment  of  rail-car  crash 


Given  the  above  assumptions,  and  assuming  a passenger  relative 


velocity  of  V at  impact,  and  a deceleration  distance  of  d,  the 
passenger  deceleration  level  (in  g-units)  is 


r 


ap  2dg  * 


The  severity  index  is: 


a 


2.5 

L 

P 
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• or 


SI 


V 

2.5  r 1.5  v 

a = a — 

P apg  p g 
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3.  TRAIN-COLLISION  PROGRAM 


3 . 1 INPUT 


Refer  to  Figure 

3-1. 

Card  1 

Format  20A4 

TITLE(I),  1 = 1,  20 

Run  title.  Up  to  80  characters  of 
alphanumeric  information  describing 
the  run 

Card  2 

F or  mat  4F  10.0, 

15 

TO 

initial-simulation  time 

TF 

final-simulation  time 

DT 

integration-time  increment 

DTP 

output -print  interval 

MODE 

MODE  = 0 assumes  series  springs 
between  cars 

MODE  = 1 assumes  inertial  masse 
acting  between  cars 

Card  3 Format  15,  5X,  5F10.0 


NCARS 


WEIGHT 

VW 

VTRN 

ACO 

APO 


number  of  cars  in  the  train  (if 
MODE  = 1,  total  number  of 
masses) 

weight  of  each  car 
barrier  velocity- 
train  initial  velocity 
train-braking  deceleration  level 
passenger  deceleration 


sec 

sec 

sec 

sec 


lb 

mph 

mph 

g 

g 
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INPUT  CODING  FORM 


cravD 


•M-  to  so  ■<  r- 

vO 


<f  00  < 

r*  ro 


CM 
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Figure  3-1  Train-Collision  Program  Input,  Sheet  1 of 


Figure  3-1  Train-Collision  Program  Input,  Sheet 


Card  4 Format  6F10.0 


SPAC(I),  I = 1,6 


passenger-spacing  distances 

before  impact  ft 


Card  5 Format  6F10.0 
CR(I ) , I = 1,5 


passenger-deceleration  distances 
after  impact  in. 


Card 

6 Format 

8G10.0 

FCO (I , 1) , 

I = 1,8 

Card 

6A  Format 

2G10.0 

FCD (I, 1) , 

I = 9,10 

Card 

7 Format 

8G10.0 

FCO (1,1)  , 

I = 1,8 

Card 

7A  Format 

2G10.0 

FCO  (1 , 2)  , 

I = 9,10 

Card 

8 Format 

8G10.0 

DEL  (1,1), 

I = 1,8 

Card 

8A  Format 

2G10.0 

DEL (1,1) , 

I = 9,10 

force  table  for  first  force- 
deflection  characteristic  lb 


force  table  for  second  force- 
deflection  characteristic  lb 


deflection  table  for  first  force- 
deflection  characteristic  ft 
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Card  9 Format  8G10.0 


DEL (1,2),  I = 1,8 

Card  9A  Format  2G10.0 

DEL (1,2),  I = 9,10 

Card  10  Format  2G10.0 
URATE(l) 

URATE (2) 

Card  11  Format  2014 

IF (I),  I = 1,  NCARS 


Card  12  Format  1G10.0  (to  be  read 
WFACTR 


deflection  table  for  second 
force-deflection  characteristic 


unloading  rates  for  first  and 
second  force-deflection 
characteristics 

force-deflection  characteristic 
number  (1  or  2)  to  be  used  for 
spring  I.  Note  that  spring  I 
acts  between  masses  1-1  and  I 

if  MODE  = 1) 

multiplier  used  to  reduce  the 
weight  of  masses  acting  between 
cars  if  MODE  = 1 
0 WFACTR  1.0. 


lb/ft 
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3.2  OUTPUT 


Output  from  the  Train-Collision  Program  includes  a listing 
of  program  inputs  as  shown  in  Figure  3-2,  and  a time  history  of  dynamic 
output  as  shown  in  Figure  3-3.  The  dynamic  output  consists  of: 

a.  The  simulated  time,  in  sec 

b.  The  car  (or  spring)  number  (note  that  car  number  0 
indicates  the  barrier  for  which  only  the  position  and 
velocity  are  output)  and  for  MODE  = 1,  even-numbered 
cars  represent  inertial  masses 

c.  XC,  the  displacement  of  each  car  relative  to  its 
position  at  impact,  in  ft 

d.  DEFL,  the  deflection  of  the  indicated  spring. 

Spring  I acts  between  masses  I- 1 and  I.  In  the  case 

of  MODE  = 0,  the  deflection  of  spring  I is 

1/2  (Xc  . - Xc  . . ) . In  the  case  of  MODE  = 1 , the 
i '-l-l  ’ 

deflection  of  spring  I is  (XQ^  - XC;j_^) 

e.  VC,  the  velocity  of  the  indicated  mass  in  ft/sec 

f.  AC,  the  acceleration  of  the  indicated  mass,  in  g 

g.  XPC,  the  relative  displacement  of  the  passenger  with 
respect  to  the  car,  for  the  indicated  car,  in  ft 

h.  VPC,  the  relative  velocity  of  the  passenger  with  respect 
to  the  car  for  the  indicated  car,  in  ft/sec 

i.  F,  the  force  produced  by  the  indicated  spring,  in  lb. 
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Figure  3-2  Sample  Listing  of  Program  Inputs 
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Figure  3-3  Sample  Program  Output 
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3.3  SAMPLE  RUNS 


Four  sample  runs  using  the  Train-Collision  Program  are  pre- 
sented in  this  section.  The  first  is  a simple  linear  spring-mass  system  used 
to  check  the  accuracy  of  the  integrator.  As  shown  in  the  program  inputs 
in  Figure  3-4,  the  mass  weighs  10  lb  and  the  spring  has  a rate  of  1226  lb /ft. 
The  exact  solution  to  this  system  yields  a natural  frequency  of  10  cps,  and 
with  an  initial  velocity  of  30  mph,  the  maximum  displacement  is  0.7003  ft, 
and  maximum  acceleration  is  -85.85  g.  The  results  from  the  program  indicate 
a maximum  displacement  of  0.70  ft  (Figure  3-5)  and  maximum  acceleration  of 
-85.85  g. 


The  second  and  third  sample  runs  compare  program  results  under 
different  modes  of  operation.  Inputs  for  the  second  run  are  shown  in  Figure 
3-6.  Twomass.es,  each  of  which  weighs  10  lb,  with  linear  springs,  of  rates 
1226  lb/ft  impact  a stationary  barrier  at  30  mph.  Displacements  of  the 
masses  are  shown  in  Figure  3-7. 

The  third  sample  run  duplicates  the  conditions  of  the  second 
(Figure  3-8)  except  that  program  mode  1 is  used.  The  two  10-lb  masses  are 
separated  by  a mass  of  0.1  lb  (WFACTR  = 0.01).  Results  of  this  run  are 
shown  in  Figure  3-9.  Note  that  the  displacements  of  the  two  10-lb  masses 
are  nearly  identical  with  the  second  run.  Although  this  mode  of  operation 
may  be  useful  in  determining  deflections  of  individual  cars  with  different 
force-deflection  characteristics,  care  must  be  exercised  in  the  choice  of 
integration  step  size  to  insure  that  stability  in  the  integration  of  the  small 
mass  dynamics  is  maintained. 

The  fourth  sample  run  is  typical  of  the  inputs  required  for  a 
four-car  train.  The  inputs  used  are  shown  in  Figure  3-10  and  the  deflections 
of  the  forward  ends  of  each  car  are  shown  in  Figure  3-11. 
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Figure  3-4  Inputs  For  Sample  Run 
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Figure  3-5  Mass  Displacement  - Sample  Run  1 
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Figure  3-6  Inputs  For  Sample  Run  2 
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Figure  3-7  Car  Displacements  - Sample  Run  2 
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TRAIN  COLLISION  MODEL 
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Figure  3-8  Inputs  For  Sample  Run 
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Figure  3-9  Mass  Displacement — Sample  Run  3 
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Figure  3-10  Inputs  For  Sample  Run 
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Figure  3-11  Car  Deflections  - Sample  Run  4 
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3.4  PROGRAM  GLOSSARY 


Program  Symbol 

Definition 

AC 

car  or  mass  acceleration 

ACO 

car-braking  deceleration 

AIMP 

passenger  impact  deceleration 

AP 

passenger  acceleration 

APO 

passenger-braking  deceleration 

CR 

passenger  crush  distance 

DEF 

car  deflection 

DEFV 

car-deflection  velocity 

DEL 

deflection  table 

DM 

deflection  start  of  unloading 

DP 

previous  value  of  DEF  (at  T-DT) 

DR 

zero  force  intercept  of  reloading  curve 

DT 

integration  increment 

DTP 

printout  interval 

DVP 

previous  value  of  DEFV  (at  T-DT) 

FC 

intermass  force 

FCO 

force  table 

FM 

force  level  at  start  of  unloading 

G 

acceleration  of  gravity,  32.  2 ft/ sec^ 

IF 

force-deflection  number  indicator 

IMP 

passenger-impact  indicator 

MODE 

indicator  for  simulation  mode  of 
operation 

3-18 


Program  Symbol 

Definition 

N,  NCARS 

number  of  cars  or  masses 

SI 

passenger  severity  index  table 

SPAC 

passenger  spacing  before  impact 

T 

simulated  time 

TF 

end  time 

TITLE 

run  title 

TO 

initial  time 

URATE 

unloading  rate 

VC 

car  or  mass  velocity 

VIMP 

passenger  relative  velocity  at  impact 

VP 

passenger  velocity 

VTRN 

train  initial  velocity 

VW 

barrier  velocity 

WEIGHT 

car  weight 

WFACTR 

factor  for  reducing  weight  of  inter - 
car  inertial  masses 

WGT 

car  or  mass  weight 

XC 

car  position 

XP 

passenger  position 

XW 

barrier  position 
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C ***  RAIL  CAR  COLLISION  MODEL  ***** 

C **  ** 

COMMON  IMP( 20) » WGT ( 20) »FC(21 ) » XC ( 20) »VC (20) , AC ( 20 ) , XP (2 0 ) , VP( 20 ) , 

1 AP(20) ,VIMP< 20,6 ) ,DVP(20) , DP ( 20 ) ,DM ( 20 ) , FM ( 20 ) , URATE (2) , 

2 DR (20) , D EL( 10»2) ,FC0( 10 ,2 ) , I F ( 20 ) ,T I TLE ( 2 0 ) ,T0,NCARS »TF , 

3 WEIGHT  * DT  »VW,DTP»VTRN»M0DE»AC0,AP0»SPAC(6)  ,CR( 5) tWFACTR, 

**  XW  ,T,G,N,AIMP(20,6,5),SI(20,6,5) 

C**  ** 

C CLEAR  COMMON  STORAGE  TO  ZERO 

C **  ** 

1 CALL  CLEAR(  IMP(1  ), SI (20,6,5)  ) 

o**  ** 

C INITIALIZE  OCCUPANT  CONTACT  INDICATOR 

C **  ** 

4  DO  2 1=1,20 

2 IMP ( I ) = 1 
C **** 

C OBTAIN  INPUT  AND  INITIALIZE 

C **  ** 

3 CALL  INPUT 
TPRNT  = TO+DTP 
G = 32.2 

FPS  = 5280. /36„0. 

VWALL  = VW  * FP  S 
VTRN  = VTRN*FPS 
N = NCARS 
N 1 = 1 
DO  31  I = 1 , N 
31  WGT { I ) = WEIGHT 
C **  ** 

C M3DE  = I I ND I CATES  THAT  EVERY  OTHER  CAR  IS  AN  INERTIAL  MASS 

C OF  WEIGHT  WFACTR+WGT  BETWEEN  CARS  SO  THAT  DIFFERENT 

C FORCE-DEFLECTION  CHARACTERISTICS  CAN  BE  USED  ON 

C ADJACENT  CARS 

C **  ** 

c **  ** 

C INITI ALIZ ATI LN 

C**** 

I F ( MODE .NE.l)  GO  TO  28 


DO  29 

I =2,N  ,2 

29 

WGT ( I ) 

= WG  T ( I ) *WFACTR 

28 

XW  = 0 

.0 

VW  = VWALL 

T = TO 
DO  32 

I =1  ,N 

FC(  I ) 

= 0.0 

XC  ( I ) 

= 0.0 

VC  ( I ) 

= VTRN 

AC  ( I ) 

= ACQ*G 

XP  ( I ) 

= 0.0 

VP  ( I ) 

= VTRN 

32 

AP(  I ) 

= 0.0 

CALL  OUTPUT(l) 

BEGIN  INTEGRATION  LOOP.  INTEGRATION  BY  MODIFIED  EULER  METHOD 
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o r o o o o o o 


VW*T 


10  xw  = 

C***« 

C EVALUATE  ACCELERATION  AT  T 

C**** 

CALL  FORCE 
DO  707  1=1, N 
ACT  = 0.0 

I F ( ABS( VC( I ) ) .GT.0.01)  ACT  = ACO+VC ( I)/ABS ( VC ( I ) J 
XDD  = ( FC ( I +1  ) — F C ( I ) )*G/WGT ( I ) +ACT 

C**** 

PROJECT  POSITION  AND  VELOCITY  AT  T+DT 

** 

XD  = VC(I)+XDD*DT 

X = XC( I)+VC(I)*DT+.5*XDD*DT*DT 
XC  ( I ) = X 

VC (I)  = XD 
707  AC ( I ) = XDD 
T = T+DT 
XW  = VW*T 

**** 

EVALUATE  ACCELERATION  AT  T+DT 

**** 

CALL  FORCE 
DO  706  1=1, N 
ACT  = 0.0 

IF(ABS(VC(I )). GT.0.01)  ACT  = ACO*VC ( I ) / ABS ( VC ( I ) ) 

XDD  = (FClI+1)  — F C ( I 1 )*G/WGT ( I )+ACT 

**** 

REVISE  POSITION  AND  VELOCITY  AT  T+DT 

** 

XD1  = V C ( I ) —AC ( I )*DT 

XD  = XD  l + .5*( AC( I) +XDD ) *DT 

XI  = XC( I )-XDl*DT-.5*AC( I )*DT*DT 

X = Xl+XDl*DT+.25*( AC( I )+XDD)*DT*DT 
XC  l I ) = X 
VC  ( I ) = XD 
706  AC ( I ) = XDD 
N 1 = 1 

IF ( MODE .EQ. 1)  N1  = 2 

C***+ 

C INTEGRATE  POSITION  AND  VELOCITY  OF  PASSENGERS  IN  EACH  CAR 

C **  ** 

DO  40  I=1»N,N1 
VDIF  = VC(I ) —VP ( 1) 

AP( I ) = 0.0 

IFIVDIF.NE.C.O)  AP(I)  = APO*VDIF/AES( VD IF  ) 

VP(I)  = VP( 1)+ AP«I )*DT 

XP(I)  = XP ( I ) +VP (I )*DT+.5*AP (I )*DT*DT 

C**** 

C CALCULATE  RELATIVE  POSITION  AND  VELOCITY  OF  PASSENGER 

C 

XPR  = XP(I)-XC (I) 

VPR  = V P ( I ) -VC  (I  ) 

IS  = IMP ( I ) 
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IF ( IS.GE.7)  GO  TO  40 

♦♦ 

IF  PASSENGER  HAS  MOVED  INPUT  DISTANCES  OF  SPAC  RELATIVE  TO 
CAR,  LOAD  RELATIVE  VELOCITY  INTO  VIMP  FOR  SUBSEQUENT  SEVERITY 
INDEX  CALCULATIONS 
SPAC  MUST  BE  IN  INCREASING  ORDER 

**** 

IF ( XPR. LT.SPAC ( IS) ) GO  TO  40 
IMP ( I ) = IS+1 
VIMP(I,IS)  = VPR 
*0  CONTINUE 

IF(T.GT.TF)  GO  TO  900 

C**  ** 

C OUTPUT  IF  PRINT  INTERVAL  HAS  BEEN  REACHED 

C **** 

IF ( TPRNT.GT ,T+ .1*DT ) GO  TO  10 
CALL  OUTPUT (2  ) 

TPRNT  = TPRNT+  DTP 
GO  TO  10 

C **  ** 

C AT  END  OF  RUN  , CALCULATE  SEVERITY  INDEX 

900  CALL  SIC 
GO  TO  4 
END 
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SUBROUTINE  OUTPUT(IND) 

C**  ** 

C PRINTS  DYNAMIC  OUTPUT  FROM  TRAIN  COLLISION  MODEL 

C**** 

COMMON  IMP( 20) ,WGT(20) ,FC(21 ),XC(20),VC (20) ,AC(20) ,XP <20 ) , VPI 20 ) , 

1 API  2 0) ,VIMP(20,6),DVP(20) ,DP 1 20 ) ,DM I 20 ) ,FM 1 20 ) »URATE!2) , 

2 DR (20),  DELI  10,2) ,FCO(10,2 ) ,IF ( 20 ) ,T ITLE I 20 ) , TO ,NC ARS *TF , 

3 WEIGHT, DT,VW, DTP, VTRN, MODE, ACO,APO,SPAC (6 ) ,CR(5) , WF ACTR , 

4 XW  ,T  ,G,N,AIMP(20,6,5),SI(20,6,5) 

C**** 

C FOR  IND=1 , CALCULATE  IPRNT  FOR  OUTPUT  PAGING 

C**  *# 

IFIIND.GT.l ) GO  TO  10 
IPRNT  = 53/ INC  ARS+2 ) 

5  LL  = 0 

C OUTPUT  PAGE  HEADINGS 

C**** 

WR I TE I 6 ,2000 ) 

2000  FORMAT! 1H1, IX , 4H TIM E ,6X , 3HC AR , 5X , 2HXC , 8 X ,4H DEFL , 6X , 2H  VC  , 

1 8X,2HAC,8X,3HXPC,7X,3HVPC,7X,1HF  / 

2 2X,3HSEC,7X,3HNG. , 5X , 3HFT. , 7X, 3HFT . ,7X ,6HFT/SEC , 4X , 3HG 'S , 

3 7X  ,3HFT.  , 7X,6HFT/SEC,4X,3HLB.  //  ) 

10  LL  = LL+1 

C*$** 

C LL  IS  A COUNTER  FOR  PAGING 

C*£** 

IFILL.GT.  IPRNT ) GO  TO  5 
I = 0 

C**** 

C OUTPUT  FOR  BARRIER 

WR ITE (6 ,200 1 ) T , I , XW, VW 

2001  FORMAT! 1H0 , F8 . 3 , 2X , I 2 , 2 X , F8  .2 , 12 X , F8 .2 ) 

FCTR  = .5 

IF ( MODE  .EQ.  1 ) FCTR  = 1.0 
DO  34  1=  1 , N 
FI  = FCII) 

XI  = XCII) 

VI  = VCII) 

A I = AC(I)/G 
XPC  = XP ( I ) -XI 
VPC  = VPID-VI 
IFI  I.GT  .1  ) GO  TO  32 
DI  = XI-XW 
GO  TO  31 

32  DI  = (Xl-XC (1-1 ) )*FCTR 
C **** 

C OUTPUT  DYNAMIC  DATA  FOR  MASSES  AND  SPRINGS 

31  WRI TE I 6 ,2002 ) I , XI , D I , V I , A I , XPC , VPC, F I 

2002  FORMAT l 1H  , 10X, I 2, 6 ( 2X , F8 .2 ) , 2X , E 10 .5 ) 

34  CONTINUE 

RETURN 

END 
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SUBROUTINE  FORCE 


C CALCULATES  FORCE  ACTING  BETWEEN  CARS  AS  A FUNCTION  OF  DISTANCE 

C**** 

COMMON  IMP ( 20) , WGT < 2 0 ) , FC ( 2 I ) , XC ( 20 ) , VC { 20 ) ,AC ( 2 0 ) ,XP (2 0 ) , VP ( 20 ) , 

1 AP(20) , VIMP(20,6 ) ,DVP<20) ,DP<20) ,DM(20) ,FM<20) ,URATE(2) , 

2 DR(20),DEL(10,2) ,FCO( 10,2)  ,IF( 20 ) , TITLE ( 2C ) ,TC ,NC ARS ,TF , 

3 WEIGHT,DT,VW,DTP,VTRN,M0DE,AC0,AP0,SPAC(6) ,CR ( 5 ) , WF ACTR , 

A XW  ,T  »G,N,AIMP(20,6,5),SI{20,6,5) 

DO  50  1=1, N 

C**** 

C NF  INDICATES  THE  APPROPRIATE  FORCE-DEFLECTION  CHARACTERISTIC 

C FOR  SPRING  I 

C **  ** 

NF  = IF  ( I ) 

IF(I.NE.l)  GO  70  51 
DEF  = XC(1)-XW 
DEFV  = VC( 1 ) —V  W 
GO  TO  52 

51  FACTOR  = l.C 

C IF  MODE  = 0 , FORCE  DEFLECTION  CHARACTERISTICS  BETWEEN  CARS 

C ARE  SERIES  COMBINATIONS  OF  THE  INDIVIDUAL  CAR  PROPERTIES. 

C THE  DEFLECTION  OF  EACH  CAR  IS  ASSUMED  TO  BE  .5*THE  TOTAL. 

C **  ** 

IF ( MODE • EQ. C ) -ACTOR  = .5 
DEF  = { XC ( I )-XC  (1-1 ) )* FACTOR 
DEFV  = VC ( T )-VC(I-l  ) 

52  IFOEF.GT.O.)  oO  TO  10 
F = 0.0 

GO  TO  300 

10  IF ( DEFV  .GT • G.  ) GO  TO  ICO 
200  IF ( DVP ( I ) . L T • C . ) GO  TO  22C 
C**  ** 

C SET  UNLOADING  CURVE 

C**»* 

DM ( I ) = DEF 
FM  ( I ) = FC ( 1 ) 

C CALCULATE  UNLOADING  FORCE 

C 

220  F = FM ( I ) — URAT n ( NF ) * ( DM ( I ) — DEF ) 

IF(F.GE.O.O)  GO  TO  300 
230  F = 0.0 

GO  TO  300 

ICO  IF ( DVP(I) .LE.Q.O)  GO  TO  150 
110  I F ( DR ( I ) .EQ.O.C  ) GO  TO  18C 
GO  TO  1 t>C 

C **  ** 

C SET  RELOADING  CURVE 

C DR  IS  THE  ZERO  FORCE  INTERCEPT  FOR  RELOADING 

C ***.$ 

150  DR ( I ) = DP(  I)-FC(I)/URATE(NF) 

C +*  ** 

C CALCULATE  FCRCE  FOR  RELOADING  ALONG  UNLOADING  CURVE 
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C**** 

160  FP  = URATE(NF)*(DEF-DR(I) ) 

IFJFP.LT. 0.0)  FP  = 0.0 

c**  ** 

INTERPOLATE  FORCE  ON  LOADING  CURVE 

**** 

180  DO  185  K=2  » 10 

IF(DEF.GT.DEL(K,NF> ) GO  TO  185 
X 1 = K— 1 

F = FC0 (K1»NF >+<FC0( K,NF )-FC0 ( K1 , NF ) )*( DEF-DEL <K 1 ,NF ) ) 

1 /( DEL( K,NF)-DEL(K1,NF) ) 

GO  TO  186 
165  CONTINUE 

F = FCO (9,NF) +( FCOJ 1 0, N F ) -FCO < 9 , NF ) )* < D EF-DEL ( 9 , NF ) )/ 

1 (DEL(10,NF)-DEL(9,NF) ) 

186  CONTINUE 

IF(DRd).EQ.O.O)  GO  TO  300 

**  ** 

IF  RELOADING,  MAKE  SURE  FORCE  LEVEL  ON  RELOADING  CURVE  DOES  NOT 
EXCEED  FORCE  ON  ORIGINAL  LOADING  TABLE 

**  ** 

IF(F.GT.FP)  GO  TO  19G 
DR ( I ) = 0.0 
GO  TO  300 
190  F = FP 
300  FC(I)  = F 
** 

SET  PREVIOS  VALUES  FOR  NEXT  TIME 

**** 

DP  ( I ) = DEF 
DVP(I)  = DEFV 
50  CONTINUE 

FCJN+l)  = 0.0 

RETURN 

END 
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SUBROUTINE  SIC 


L **  ** 

C CALCULATES  SEVERITY  INDEX  MATRIX  FOR  EACH  CAR  FOR  PASSENGER 

C DISTANCES  OF  SPAC  BEFORE  IMPACT  AND  STOPPING  DISTANCES 

C OF  CR.  NOTE  THAT  THE  RELATIVE  VELOCITY  OF  THE  PASSENGER 

C W/TG  THE  CAR  AT  SPACING  DISTANCES  SPAC  ARE  VIMP  AS  CALCULATED 

C IN  MAIN 

C*** 

COMMON  IMP( 20) , WGT(2C) , FC<21 ) ,XC (20) »VC (20) ,AC(20 ) , XP (20 ) ,VP(2C ) , 

1 AP(2G) ,VIMP< 20»6) ,DVP(20) ,DP(20) »DM(20) ,FM(20) ,URATE(2), 

2 DR (20) , DEL ( 10,2 ) ,FCO< 10,2  ) , I F ( 20 ) , TITLE ( 20 ) , TO ,NC AR S ,TF , 

3 WEIGHT,  L>T,VW,  DTP,  VTRN,  MODE  , ACO  , A PO , S PAC  ( 6 ) ,CR(  5)  , WF  ACTR  , 

A XW  ,T  ,G,N, AIMP (20,6,5), SI<20,6 ,5) 

N 1 = 1 

IF (MODE.  tO.  1)  N 1 = 2 
DO  6 00  1 = 1, N,N  I 
DO  800  K= 1 , 6 
DO  800  M= 1 , 8 
C 

C AIMP  IS  CONSTANT  ACCELERATION  FROM  VELOCITY  VIMP  IN  DISTANCE  CR. 

C**  ** 

AIM  P ( I , K ,M ) = 6.*VIMP(I ,K)**2/(CR(M)*G) 

C SEVERITY  INDEX  IS  INTEGRAL  OF  ACCE LER ATI  ON** 2 .5  OVER  TIME 

S I ( I , K , M ) = AIMP(I,K,M)**1.5*VIMP(I,K)/G 
SCO  CONTINUE 
C **  * 

C OUTPUT  SEVERITY  INDEX  TABULATION  AS  A FUNCTION  OF  PASSENGER 

C SPACING  AND  CRUSH  DISTANCE  FOR  EACH  CAR 

C **  ** 

WR I TE ( 6 ,3000 ) 

3000  FORMAT ( 1H 1 , //  4cX, 14H SEVERITY  INDEX  ) 

DO  950  I=1»N,N1 

WR I TE ( 6 ,300 1 ) I , (SPAC(K) ,K=1,6) 

30C1  FORMAT ( 1H0, 1GX, 3HCAR ,13  / 3 X , 5HC RUSH , 22 X , 1 2HS P AC  I NG  - FT  / 

1 AX,3HIN . ,3X,oF10.2  ) 

DO  950  K= 1 , 5 

WRITE(6 ,300  2 ) CR(K),(SI(I,J,K),J  = 1,6) 

3002  FORMAT< 1H0, F7.2 ,2X, 6FI0.1  / ( IX  , F7 .2 , 2 X ,6 F 10 . 1 ) ) 

950  CONTINUE 

I F ( MODE • EQ . 1 ) WRITE ( 6,3C03) 

3003  FORMAT ( 1H0, //47H  NOTE:  M0DE=1,  THEREFORE  PASSENGERS  ASSOCIATED 

1 20HWITH  EVERY  OTHER  CAR  ) 

RETURN 

END 
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SUBROUTINE  INPUT 


C READS  TRAIN  COLLISION  MODEL  INPUT  FROM  CARDS  AND  PRINTS  INPUT  DATA 

C ** 

COMMON  IMP( 20) ,WGT(20) , FC ( 2 1 ) ,XC ( 20 ) , VC ( 20 ) , AC ( 20 ) ,XP (20 ) , VP< 20 ) , 

1 AP(2C)  ,V IMP (20,6 ) ,DVP (20) , DP (20 ) , DM ( 20 ) , FM ( 20 ) , URATE (2) , 

2 DR(2G),DEL(10,2) ,FC0( 1C  , 2 ) , I F ( 20  ) ,TI  TLE  ( 20  ) ,T0  ,NCAR  S ,TF  , 

3 WEIGHT, DT»VW»DTP,VTRN»M ODE, AC 0»APQ,  SPAC ( 6 ) ,CR( 5) ,WFACTR, 

A XW  ,T,G,N,AIMP(20,6,5),SI(20,6,5) 

C READ  INPUT 

C **  ** 

C **  ** 

C RUN  TITLE  CARD,  80  CHARACTERS 

READ(5,99)  (TITLE( I ) ,1=1,20) 

99  FORM AT ( 20A4 ) 

C**** 

C PROGRAM  CONTROL 

£♦*** 

READ (5,98)  T 0 , T F , DT , DT P , MOD E 
98  FORMAT( 4F10 .0, 15) 

READ ( 5, 97)  NC A RS  , WE  I GHT , VW , VTRN , ACO , A PO 
97  FORMAT! 15, 5X,5F10.0) 

RE A D ( 5 » 96 ) ( SPA C ( I ) , 1= 1 , 6 ) 

RE A D ( 5 , 96 ) ( C R ( 1 ) , I = 1 , 5 ) 

96  FORMATt 6F10  .0) 

**** 

EACH  F CO  AND  DEL  REQUIRES  2 CARDS  WITH  2 VALUES  ON  SECOND  CARD 
BLANK  CARD  MUST  BE  FURNISHED  IF  NO  DATA  FOR  LAST  TWO  ENTRIES. 

READ(5,95)  ( FCO (1,1), 1=1, 10) 

READ(5»95 ) ( FCO ( I , 2 ) , I = 1 , 10 ) 

RE  AD ( 5 , 95 ) ( DEL ( I , 1 ) , I = 1 , 1 0 ) 

READ( 5,95)  ( DE L ( 1 , 2 ) , I =1 , 10  ) 

RE  AD  ( 5 , 95  ) UR  A IE  ( 1 ) , UR  A T E ( 2 ) 

95  FORMAT! 8G10 .0) 

** 

IF  INDICATES  WHICH  FCQ-DEL  TABLE  USED  FOR  EACH  SPRING.  MUST 
BE  1 OR  2 

**  ** 

RE AD( 5, 94 ) ( I F ( I ) , I = 1 , NC ARS  ) 

94  FORMAT ( 2014  ) 

C WFACTR  READ  ONLY  IF  MODE  = 1 

IF ( MODE . EQ. 1 ) READ  ( 5,95 ) WFACTR 

O*** 

C PRINT  INPUT 

C **  ** 

WRITE (6 ,1000)  ( TITLE! I i ,1=1,20) 

1000  FORMAT( 1H1, 40X, 2 1HTRAIN  COLLISION  MODEL  / 20X,2GA4) 

WRITE (6, 1001)  T G,NC ARS, TF, WE IGHT,DT,VW, DTP, VTRN, MODE, ACO, APO 
WRITE (6, 1002)  ( FCO( I ,1)  , DEL (1,1) , FCO( 1 , 2 ) , DEL ( 1,2) ,1  = 1,10) 

WR I TE ( 6 , 100  3 ) URATEt 1) ,URATE(2) 
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WRITE (6 , 1004)  ( 1,1  = 1, NCARS) 

WRITE(6,1005)  ( I F ( I ) ,1=1, NCARS) 

WRITE (6 ,1006)  ( SPAC ( I) ,CR( I ) ,1=1, 5) ,SPAC(6 ) 

I F ( MODE • EQ • 1 ) WRITE ( 6, 1007)  WFACTR 
10C7  FORMAT ( 1H0,  //  9H  WFACTR  =,F10.4) 

1001  FORMATt lH0,25Xt^HTC  =,F8.4,4H  S EC , 8X , 13HN0 • OF  CARS  =,14,  / 

1 26X , 4HTF  =,F8.4,4H  S EC ,9X , 1 1HCARWE IG HT  =,F12.1,4H  LBS  / 

2 26X,hHDT  = , F 8 .5 , 4H  S EC , 3X , 18HB ARRI ER  VELOCITY  =,F12.2,4H  MPH/ 

3 20X,  10HDT  PRINT  =,F8.4,4H  SEC  » 5X , 16H TRAIN  VELOCITY  =,F12.2, 

4 4H  MPH  / 

5 24X , 6HM0D  E = , 1 4, 1 5 X , 14HTRA IN  DECEL.  =,F12.2,4H  G'S  / 

6 45X, 13HPASSENGER  DECEL.  =,F12.3,4H  G'S  ) 

1002  FORMAT ( 1HC » //  1QX,38HTRAIN  FORCE-DEFLECTION  CHARACTERISTICS  / 

1 13X,6HFC0( 1) ,10X,6HDEL(1),10X,6HFCO(2),10X,6HDEL(2)  / 

2 14X,3HLBS,13X,3HFT. , 13X, 3HLBS , 13X, 3HFT.  // 

3 ( 10X,F10.1,3X,F7.2,7X,F10.1,8X,F7.2)  ) 

1003  FORMAT( 1H0, 10X, 26HUNL0ADING  RATES:  UR AT E ( 1 ) = , F 14 . 2 , 5 H LB / FT , 

2 10X,9HURATE( 2 )=, F14.2, 5HLB/FT  /) 

1004  FORMAT ( 1H0,  //  5X,7HCAR  NO. ,23X, 1516  ) 

1005  FORMATt 1H  , 4X , 3 5HF0RCE-DE F L E CT ION  CHARACTERISTIC  NO.  ,1515  ) 

lOGo  FORMAT ( 1H0 , //  10X, 17H PASSENGER  S P AC I NG , 1 3 X , 16HP A S SE N GE R DECEL.  / 

1 10X, 13HDISTANCE  - FT , 17X , 13HDI STANCE  - IN  / 

2 ( 13X ,F8. 2,22X,F8 .2)  ) 

RETURN 

END 
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SUBROUTINE  CLEAR ( A » B ) 

CLEARS  (SETS  TO  ZERO)  A BLOCK  OF  STORAGE  IDENTIFIED  BY 
THE  ADDRESSES  OF  THE  TWO  ARGUMENTS 

CALL  CLEAR ( PtQ ) WILL  CAUSE  ALL  BYTES  TO  BE  SET  TO 
ZERO  FROM  ADDRESS  P THROUGH  THE  FULL-WORD  ADDRESS  AT  Q 
DIMENSION  A ( 1 ) ,B(1) 

3(1)  = 1.0 
1 = 0 

10  IF ( B( 1) .EQ.O.O)  RETURN 
I = I + l 
A ( I ) = 0.0 

GO  TO  10 
END 
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3.6  PROGRAM  FLOWCHARTS 
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AUTOFLOW  CHART  SET  ~ — TRAIN  COLLISION  MODEL 


CHART  TITLE  - PROCEDURES 


•••  «AIL  CAR  collision  mooel 


I 

CLEAR  COMMON  STORAGE  TO  ZERO 
I 01 


I 1 I 

I 6 | CLEAR 

I.  I ( IMP  ( 1 1 * SX  «*2C  • 

I 0 I 6,  i> ) ) 

I 1 I 


I 

INITIALIZE  OCCUPANT 
CONTACT  I NDICATO  * 

I 

04.12 >1 

*♦  I NOT  * 02 

* BEGIN  DO  LOO p * 

* 2 I = 1,  20  * 


IMPI  1 ) = 1 


* END  OF  DO 

• LOOP? 


I 

I 

OBTAIN  INPUT  AND 
INITIALIZE 

I 

I 

3 I 


II  I 

13  | 

i.  r 
10  | 
II  I 


TPRNT  = TO  ♦ OTP 
G = 32.2 

FPS  = 5280. /36C0. 
VWALL  = VW*FPS 


VTRN  = VTR  N*F  PS 
N = NCARS 


I NOTE  00 

* BEGIN  DO  LOOP  * 

♦ 31  I = It  N * 

I 

311  09 

I WGTIII  = WEIGHT  | 

I 


>* 

I 


NO  * 

* ENO  OF  DO 

• LOOP? 


->* 


MOD  E = 1 INDICATES  THAT 
EVERY  OTHER  CAR  IS  AN 
INERTIAL  MASS 
OF  WEIGHT  WF ACTR *WGT 
BETWEEN  CARS  SO  THAT 
DIFFERENT 
F0RCE-0EFLECT10N 
CHARACTERISTICS  CAN 
BE  USED  ON 
ADJACENT  CARS 


INITIAL 1ZA 


* MODE  .NE.  1 


* BEGIN  DO  LOOP 

* 29  I x 2,  N,  2 

I 

C2.14 >| 

29  | 

I WGTIII  = 

I WGT ( I I *WF ACT R 

I 


END  OF  DO 
LOOP? 


NO 

*-♦ 

I 


I NOTE 

******••• 

* BEGIN  DO  LOOP 

• 32  I * 1.  N 


FCIII  = 0.0 
XC(  I)  = 0.0 
vein  * VTRN 
ACI  I I « AC 0 *G 
XPI 1)  = 0.0 
VPI I)  = VTRN 


32  I 

API  I) 


END  OF  DO 
LOOP? 


* • FALSE 

• AB  SI  VC  (I)  I * 

* .GT.  C.Ol  • 


I ACT  = | 

I ACO*VC(I>  I 

I /ABS I VC  I 1 1 I | 

* — * 

I 

I 27 

| XDD  = IFCII  ♦ | 

I II  - I 

I FCII I )*G/WGT(1 ) ♦ | 
I ACT  | 

I 


PROJECT  POSITION  AND 
VELOCITY  AT  T »DT 

I 


XD  = VC  1 1 I ♦ 
XDD*DT 

X = XC  II  I ♦ 
VC  I I l*DT  ♦ 
•5*XDD*DT*DT 

XCII)  = X 

vein  = xd 


707  | 

ACtl)  = XDD 


END  OF  DO 
LOOP? 


I T = T 

I 

I XW  = 


EVALUATE  ACCELERATION 
AT  T+DT 

I 

I 
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AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLfc  - PROCEDURES 


02. 

32 >• 

♦- 

— 

>♦ 

1 

01  1 

| 

* 



— * | 

* 10 

| | 

H | 

• • 

18  | 

FORCE 

H | 

• * 

1 . 1 

H | 

NO  * * 

10  1 

H | 

* END  OF  DO  * 

1 1 1 

H | 

1 

* LOOP?  * 

* 

— 

* 1 

i 

* • 

1 

1 

* * 

1 

i 

• 

I NOTE 

02  1 

1 

1 YES 

* 

BEGIN  DO  LOOP 

• | 

1 

1 

• 

706  I = I,  N 

• | 

1 

1 

1 

1 

1 

|< 

1 

1 

1 

1 11 

1 

03  1 

• * 

*— — 

— 

— * 1 

1 N 1 = 1 1 

1 

ACT  = 0.0 

1 1 

• * 

* 

1 

— * | 

| 

1 

1 

* 12 

• 04 

• * 

• • 

• • 

* • 

FALSE  • * 

ALSE 

• • 

* MODE  .EQ.  1 * 

* 

ABS( VC( I ) ) 

• | 

1 

• * 

• .GT.  0.01  • 

1 

• • 

* • 

1 

• • 

• * 

1 

* 

• 

1 

1 TRUE 

ITRUE 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 13 

| 

05  1 

1 

« • 

• 

— — — 

— * 1 

1 

1 N 1 = 2 1 

| 

ACT  = 

1 1 

1 

• * 

1 

AC0*VC(  I ) 

1 1 

1 

| 

1 

/ A BS ( VC ( I ) ) 

1 1 

> j 

| 

INTEGRATE  POSITION 

— 

>| 

AND  VELOCITY  OF 

1 

06  | 

PASSENGERS  IN  EACH 

• 

CAR 

1 

XDD  = (FC(I  ♦ 

| | 

1 

1 

1)  - 

1 1 

1 

I FC  ( I ) )*G/WGT ( I ) 
I ACT 


REVISE  POSITION  AND 
VELOCITY  AT  T ♦DT 

I 


XD1  = VC (II  - 
AC( I ) *DT 

XD  = XDI  ♦ 

• 5* ( AC  I I ) ♦ 
XDD  J *DT 

XI  = XCII)  - 
XDI*DT  - 
. 5*  AC ( I l*DT*OT 


X = XI  ♦ XD I*DT  ♦ 

• 2 5* ( AC ( I ) ♦ 

XDD) *DT*DT 

XCII)  = X 

veil)  = XD 

I 

I 

706  I 0 

| ACID  = XDD 


I NOTE 

• •***•** 

BEGIN  DO  LOOP 
40  I = I, 

I 

<*.06 >1 


VDIF  = VC<1> 
VP(  I ) 


N1 


* * FALSE 

• VDIF  .NE.  0.0  *-♦ 

* * I 


I 

I 17 

| A P C I I = I 

| APO*VDIF/ABS  I 

| (VDIF)  I 
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C3.I6* — >* 

I 01 

| VP(  I)  = VP  I I)  ♦ I 
APII|*OT  I 

XP(I)  * XP(I)  ♦ I 
VP(I»*DT  ♦ | 

•5*AP( I»*DT*DT  | 


I 

I 

CALCULATE  RELATIVE 
POSITION  ANO  VELOCITY 
L'F  PASSENGER 

** 

I 

I 

I 02 

| XPR  = XPCI)  - I 

l xcm  I 

I I 

| VPR  = V P ( I I - | 

I vein  l 

I I 

I IS  = IHPIll  I 


-04. 03* — >* 

I 


ENO  OF  DO 
LOOP? 


NO 

I 

I 


YES  . 3 . 

. 13  . 


FALSE 


• 07 

• * 


IF  PASSENGER  HAS 
MOVED  INPUT  DISTANCES 
:F  SPAC  RELATIVE  TO 
CAR,  LOAD  RELATIVE 
VELOCITY  INTO  VIMP 
FOR  SUBSEQUENT 
SEVERI TY 

INDEX  CALCULATIONS 
»>P  AC  MUST  BE  IN 
INCREASING  ORDER 

I 

I 

I 

* C4 


* * TRUE  | 

* XPR  .LT.  * ♦ 

* SPAC  1 1 S I * A 


FALSE 


G5 


I IMP< I ) = IS  ♦ 1 I 

I I 

I VIMPU.IS)  = VPR  | 


I 

I 


* * TRUE 

* T ,GT.  TF  * ♦ 


• 

I FALSE 
I 
I 
I 

OUTPUT  IF  PRINT 
INTERVAL  HAS  BEEN 
REACHED 

I 

I 

I 

* 08 

• * 

* * 

* * TRUE 

* TPRNT  .GT.  *-♦ 

* T ♦ .1*01  • I 

* • | 

* * | 


AT  END  OF  RUN  , 
CALCULATE  SfcVER  ITY 
INDEX 


>* 

9 CO  | 11 


* — — -• 

111  H 

II  I SIC  H 

I.  I H 

10  | H 

III  H 


I 

I 

I 12 


1 FALSE 

. 2 , 

j 

. 21 

] 

10 

1 

1 

09 

* — 

1 

1 

H 

16 

| OUTPUT 

H 

1 • 

1 (2) 

H 

10 

1 

H 

1 1 

1 

H 

1 

1 

1 10 

1 

TPRNT  = TPRNT 

♦ 1 

DTP 

I 

I 


. 2.21. 

...  10 


3-34 


JV/11/  /<, 


AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - NON-PROCEDURAL  STATEMENTS 


COMMON  IMP (20) , WGT ( 2 0 ) , FC ( 2 1 ) , XC ( 20  ) ,VC ( 20  I ,AC ( 2 0) »X  P ( 20 ) , VP ( 2C > , 
AP(20). VIMP(20.6) ,DVP(20),DP(20) ,DM(20 I ,FM(20) ,URATEC2  >, 
DR120), DEL ( 10*2) #FCO( 10,2) t IF (20) .TITLE (20) ,T0 , NC AR S * T F , 

WE IGHT,DT,VW,DTP,VTRN»MODE,ACOf AP0»SPAC(6 ),CR<5 ) , WFACTR, 

XW  ,T,G,N, AIMP( 2 0.6,5 ) ,SI I 20,6,5 ) 


3-35 


09/11/74 


AUTOFLOW  CHART  SET  .MU- 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - SUBROUTINE  OUTPUT! IND) 


/ OUTPUT  / 


02.20 ♦ — >• 

PRINTS  DYNAMIC  OUTPUT 
FROM  TRAIN  COLLISION 
MODEL 


FOR  INDM*  CALCULATE 
IPRNT  FOR  OUTPUT 
PAGING 


TRUE  * 

♦ * IND  .GT.  1 


I 02 

| IPRNT  = | 

I 55/INCARS  ♦ 21  | 


I 

I 

OUTPUT  PAGE  HEADINGS 

I 


I 04 


/ / 
/ WRITE  TO  DEV  / 

/ 6 / 

/ VIA  FORMAT  / 

/ 2000  / 


I 

> | 

10  I 05 


I LL  = LL  ♦ 1 | 


I 

I 

LL  IS  A COUNTER  FOR 
PAGING 

I 

I 

I 

* 06 

* * 


* LL  .GT.  IPRNT  *- 


I FALSE 

I 

I 

I 


I 07 


I = 0 


I 


I 

I 

I 

UTPUT  FOR  BARRIER 
I 
I 


I >l  09 

I / WRITE  TO  DEV  / 

I / 6 / 

I / VIA  FORMAT  / 

I / 2001  / 

I / FROM  THE  LIST  / 

I I 

I I 

I | NOTE  10 


• LIST  « T , I,  XWt  * 

• VW  * 


I 

I 

I 11 


1 FCTR 

= .5 

1 

1 

1 

1 

• 12 

• 

• 

* 

• 

* 

* 

FALSE 

* MODE 

.EQ.  1 

• — ♦ 

• 

• 

1 

• 

• 

1 

• 

• 

1 

• 

1 

1 TRUE 
1 

1 

1 

1 

1 

1 

1 

1 

1 

13  1 

I FCTR 

* 1.0 

1 1 

1 i 

| NOTE  14 

• •*••**•*** 

* BEGIN  DO  LOOP  • 

• 34  I * It  N • ♦ 

*»**»»****•  | 

I I 

06.22 >|  I 

I 15  I 

I FI  * FCIII  I I 

I I I 

I X!  tXCIII  I I 

I I I 

I vi  s vein  i i 

I I l 

I A I 3 ACIII/G  I I 

I I I 

| XPC  = XP(I>  - XI  | | 


I 

I 

I 16 

* • 

| VPC  = VP(I)  - VI  | 

I 


* 17 


* TRUE  | 
I .GT.  1 • ♦ 

• * 

• • 

IFALSE 


>♦ 

32  | 19 

I D I » (XI  - XCII  - | 

I 1))*FC1R  | 


I 

I 

I 

OUTPUT  DYNAMIC  DATA 
FOR  MASSES  AND 
SPRINGS 

I 

>| 

31  | 20 

/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 2002  / 

/FROM  THE  LIST  / 


I 

I NOTE  2 


* LIST  = I,  XI,  DI* 

* VI,  AI,  XPC,  VPC, 

* FI 


I 

I 

I 

34  • 22 

• * 

• * 

• * NO 

• fcND  OF  DO  *-♦ 

* LOOP?  • | 

* • I 

* * | 


I DI  = XI  - XW  | 


I 

I 


YES 


23 


* EXIT  • 


3-36 


09/ll/7<» 


AUTQFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - NON-PROCEDURAL  STATEMENTS 


2000 


2001 

2002 


COMMON  IMP! 20) *WGT ( 20 ) *FC (2 1 I » XC ( 20 ) • VC ( 20 ) , AC ( 20) ,XP ( 20  ) » VP ( 20 ) « 
AP(20),VIMP(20,6) .DVP(20)tDP(20),DMI20),FM(  20 >, URATE! 2) , 

DR  I 20) , DEL CIO, 2).FC0! 10 *2) » IF (20) .TITLE (20) , T 0 , NC AR S , T F , 
WEIGHT.DT.VW.OTP, VTRN,MODE,ACO,APO.SPAC(6)*CR (5) .MFACTR, 

XW  ,T,G,N,AIMP( 20,6,5 ) ,SI( 20,6,5  I 
FORM  AT ( 1H 1 » 1X,4HTIME,6X, 3HC AR, 5 X, 2HXC . 8X , 4HDE FL . 6X ,2HVC» 
8X,2HAC*8X,3HXPC.7X, 3HVPC .7X.1HF  / 

2X  » 3HS EC » 7X  » 3HN0 • » 5X ,3HFT . #7X • 3HFT. • 7X ,6HFT /S EC • AX , 3HG*S ♦ 
7X,3HFT.,7X,6HFT/SEC,4X,3HLB.  //  ) 

FORMAT ( 1H0,F8.3,2X,I2,2X, F8.2, I2X.F8.2) 

FORMAT ( IH  , 10X.I2.6l2X.F8  J2.  ) .2X.E10.5) 


3-37 


09/11/74 


AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - SUBROUTINE  FORCE 


/ FORCE  / 


02,22* — >* 

CALCULATES  FORCE 
ACTING  BETWEEN  CARS 
AS  A FUNCTION  OF 
0 ISTANCE 


I NOT  k 01 
*********  ** 

* BEGIN  DO  LOOP  * 

* 50  I = 1,  N ♦ 

***•*••*•*• 

I 

09.17 >| 

I 

NF  INDICATES  THE 
APPROPRIATE 
FORCE-OEFLECTION 
CHARACTERISTIC 
FOR  SPRING  I 
I 
I 

I 02 


I 


NF  = IF  ( I ) 


I DEF  = XCIll  - XW  | 

I I 

I DEFV  = VCllJ  - VW  | 


I 

IF  MODE *0 , FORCE 
DEFLECTION 
CHAR ACT  ERISTICS 
BETWEEN  CARS 
ARE  SERIES 
COMBINATIONS  OF  THE 
INDIVIDUAL  CAR 
PROPERT  IES. 

THE  OEFLECTION  OF 
EACH  CAR  IS  ASSUMED 
TO  BE  • 5*THE  TOTAL. 
I 

I 


* * TRUE 

* DEFV  .GT  . 0.  *- 


* • 

IFALSE 


DVP ( I ) .LT. 


SET  UNLOAD 


I DMI I)  = DEF 
| FM  ( I > *=  FC(I  ) 


I 

I 

CALCULATE  UNLOAD ING 
FORCE 

I 

> , 

220  | 1 


I F « FM(I)  - 
| URATEINF )* (DM(I) 
I - DE  F » 


IFALSE  . 


TRUE 

*-♦ 

I 

I 

9 . 
15  . 


* * TRUE 
DVPC1I  .LE.  *-♦ 

* 0.0  ♦ | 


01  . 
150 


* * TRUE 
DR  C I ) .EU.  *-* 

* 0.0  * | 

* * I 

* * I 


FALSE 


. 9.02. 
...  160 


9.15. 
...  300 


3-38 


09/ll/7<» 

CHART  TITLE  - SUBROUTINE  FORCE 


AUTOFLOW  CHART  SET  


TRAIN  COLLISION  MODEL 


SET  RELOADING  CURVE 
DR  IS  THE  ZERO  FORCE 
INTERCEPT  FOR 
RELOADING 

08.17 >• 

150  I 01 

I DR(  I)  = DP (I)  - I 
I FCm/URATEINF)  | 


CALCULATE  FORCE  FOR 
' ELOAOING  ALONG 
UNLOADING  CURVE 

08.18 >| 

160  | 0 


I URA T E ( NF ) • ( DEF  - 
I DR  (IH 


FALSE  * 

» • FP 


I 


FP  = 0.0 


I 


I 

I 

INTERPOLATE  FORCE  ON 
LOADING  CURVE 

I 

08.18 > | 

180  | NOTE  05 

• BEGIN  DO  LOOP  • 

• 185  K = 2,  10  • 

• *••••***** 


• • TRUE 

DEF  .GT.  • 

* DEL(K,NF|  * 


END  OF  DO 
LOOP? 


I YES 

I 


I 

9 

06 


F = FC0(9,NF)  ♦ 
(FCO(lOtNF)  - 
FC0I9.NF I )• ( DE  F 
DE  L ( 9 »NF  ) ) 

/ ( DEL ( 10«NF ) - 
DE  L (9  * NF  ) ) 

I 

> | 

186  | NOTE  10 


DR  ( I ) ,E0. 

* 0.0 


IF  RELOADING*  MAKE 
SURE  FORCE  LEVEL  ON 
RELOADING  CURVE  DOES 
NOT 

EXCEED  FORCE  ON 
ORIGINAL  LOADING 
TABLE 

I 

I 

I 

* 12 


1 

07 

1 

• * 

1 

1 

50 

* 17 

* 

— 

— * 

1 

* * TRUE 

1 

1 

* 

• 

1 

K1  = K - 1 

1 

1 

• F .GT.  FP  * 



♦ 

* 

• 

1 

1 

1 

* • 

1 

* 

* NO 

1 F 

= FCO(Kl.NF)  ♦ 

1 

1 

• * 

1 

* END 

OF  CO 

1 

( F CO  ( K *NF  ) - 

1 

1 

• * 

1 

* LOOP?  * | 

1 

FC0(K1,NF I I 

1 

1 

* 

1 

• 

* 1 

1 

• ( DEF  - 

1 

1 

I FALSE 

1 

• 

• 1 

I 

DEL ( K 1 »NF  ) » 

1 

1 

1 

1 

* .... 

1 

/(DEL(K.NF)  - 

1 

1 

1 

1 

1 YE  S . 3 . 

1 

OELIKI.NFI) 

1 

— * 

1 

1 

1 

1 

1 

1 . 02  . 
1 .... 

! 

1 

1 

1 

1 13 

1 

1 

1 

i 

I DR  ( 1 1 = 0.0  1 

1 

1 

m 

1 18 

SET  PREVIOS  VALUES 
FOR  NEXT  TIME 


DP  II»  = DEF  | 

I 

DVP(I)  = DE  FV  | 


I FC(N  ♦ II  = 0.0 

* — 

I 

I 

I 19 


3-39 


J 


09/11/74 


AUTOFLOW  CHART  SET  •*' 


TRAIN  COLLISION  HODEL 


CHART  TITLE  - NON-PROCE OUR AL  STATEMENTS 


COMMON  IMP( 20) , WGT ( 20 ) »FC  12 1) , XC ( 20 ) , VC ( 2 0 ) , AC  I 2 0) ,X P{ 20 ) , VP ( 2C ) , 
API  201 ,V IMP (2 0,6 1 »DVP(20)tDP(20)tDM(20) »FM( 20 ) , UR  AT E ( 2 ) , 
DRC20>,DEL(10,2l,FC0t 10,21 , IFC20I .TITLE (20) ,T 0 ,NC AR S, T F , 
WEIGHT,DTf VW,DTP,VTRN,MOOE,ACOt APO.SPAC (6 ),CR  <51 .WFACTR, 

XW  ,T,G,N,AIMPJ20,6,5I ,SI ( 20,6,5) 


3-40 


09/11/74 


AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - SUBROUTINE  SIC 


/ SIC  / 


04.11 >* 

CALCULATES  SEVERITY 
INDEX  MATRIX  FOR  EACH 
v.AR  FOR  PASSENGER 

distances  of  spac 
pEFORE  impact  and 
STOPPING  DISTANCES 
F CR.  NOTE  THAT  THE 
RELATIVE  VELOCITY  OF 
THE  PASSENGER 
W/TO  THE  CAR  AT 
^PACING  DISTANCES 
jPAC  ARE  VIMP  AS 
CALCULATED 
IN  MAIN 


FALSE  • * 

♦ • MODE  .EQ.  1 


AlMP  IS  CONSTANT 
ACCELERATION  FROM 
VELOCITY  VIMP  IN 
DISTANCE  C 

I 

I 

I 07 

I AIMPI  I , K , M ) = | 

| 6 • *V IMP ( I , K ) I 

I **2/<CRIM)*G)  | 


I 

I 

SEVERITY  INDEX  IS 
INTEGRAL  OF 
ACCELERATION* *2.5 
OVER  TIME 

I 

I 


I SIUfKtM)  = 1 

I AI  MP  ( I *K  * M ) | 

I **1.5*VIMP( I ,R)/G  I 

I 

I 


END  OF  DO 
LOOP? 


• * NO 

• END  OF  DO 

• LOOP?  • | 

* * I 

* * I 

* .... 

I YES  . 11  . 


• 

_ 

— • 

1 

1 

• 12  1 

1 

1 

1 Nl  = 2 

1 

1 

1 

* * 1 

1 

9 50 

• 20 

• 

— — 

— • 

1 

1 

* * 1 

1 

• • 

| 

1 

1 

• * NO  | 

1 

• • 

> j 

1 

1 

• END  OF  DO  *-♦  | 

1 

NO  * 

• 

1 

NOTE 

04 

1 

1 

• LOOP?  • 1 1 

♦— 

* 

END  OF  DO 

1 

1 

* * 1 1 

• 

LOOP?  * 

• BEGIN  DO 

LOOP 

• 

1 

1 

• • II 

• • 

• BOO  I = 1, 

N,  Nl 

• 

1 

1 

• ....  1 

• • 

1 

1 

1 YES  . 11  . 1 

• 

| 

1 

1 

1 . 05  . 1 

1 YES 

11.12 >1 

1 

1 

1 ....  1 

1 

1 

NOTE 

05 

1 

1 

i i 

1 

* BEGIN  DO 

LOOP 

• 

1 

1 

1 1 

1 

* 800  K = ; 

, 6 

• 

1 

1 

i i 

1 

1 

1 

OUTPUT  SEVERITY  INDEX  | 

1 

1 

1 

1 

TABULATION  AS  A | 

* 21 

11.11 >1 

1 

1 

FUNCTION  OF  PASSENGER  1 

• * 

1 

NOTE 

06 

1 

1 

SPACING  AND  CRUSH  | 

• • 

1 

1 

DISTANCE  FOR  EACH  CAR  | 

NO  * 

• 

* BEGIN  DO 

LOOP 

• 

1 

1 

1 1 

♦— 

* 

END  OF  DO 

* 800  M = : 

, 5 

* 

1 

1 

1 1 

1 

• 

LOOP?  • 

i 

1 

1 13  ] 

1 

• * 

!<- 

— 

— 

— « 

/ / 

1 

• 

/ WRITE  TO  DEV  / 

/ 6 / 

/ VIA  FORMAT  / 

/ 3000  / 

I 

I 

I NOTE  14 
• *••••**••• 

* BEGIN  DO  LOOP  * 

* 950  I = 1,  N,  N1  * 


/ WRITE  TO  DEV 
/ 6 
/ VIA  FORMAT  / 
/ 3001  / 

/ FROM  THE  LIST  / 


| MDTE  16 

• LIST  * It  • 

• (SPAC(K).K  = 1,6)  * 

I 

I NOTE  17 

*********** 

• bEGIN  DO  LOOP  * 

• 950  K = 1,  5 * 

I 


/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 3002  / 

/ FROM  THE  LIST  / 


I NOTE  19 
»*•****•*« 

LIST  « CR  (K  ) , * 

ISI(  I ,JfK)  ,J  = * 

1,6  ) • 


* MODE  .EQ. 


/ WRITE  TO  DEV  / 


/ VIA  FORMAT  / 

/ 3003  / 


I 

I 24 


* EXIT  * 


3-41 


0<>/ll/7<* 


AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - NON-PROCEDURAL  STATEMENTS 


3000 

3001 

3002 

3003 


COMMON  IMP  I 20)  , WGT ( 20 ) ,FC 12 1 ) , XC ( 20 ) , VC ( 20 ) , AC ( 2 0 ) ,X P( 20  > , VP < 20 1 , 
AP(20),VIMP(20,6) , DVP (20),DP<  20) ,DM(20I ,FM( 20) ,URATE<2  > , 
DR(20)»DEL(10»2)»FCO(10»2)»IF(20)»TITLE(20)*TO»NCARS»TF» 
WEIGHT,DT,VW,DTP,VTRN,MODE, ACO, AP0,SPAC!6 ),CR  <51 , WF ACTR, 

XW  ,T,G,N, AIMP( 20, 6,5) ,SI < 20,6,5) 

FORM  A T ( 1H 1 , //  *0X  * 14H  SEVERITY  INDEX  ) 

FORM  A T ( 1H0  * 1 OX • 3HCAR  *13  / 3 X , 5HCRUSH , 22X , 12HSP AC ING  - FT  / 

<*X.  3HIN.  ,3X,6F10.2  1 

FORMAT! 1H0»F7.2»2X,6F1C.1  / I IX, F7 .2 * 2X , 6F 10 . 1 ) > 

FORMAT | lHO* //4  7H  NOTE:  MODE-1,  THEREFORE  PASSENGERS  ASSOCIATED 
20HWITH  EVERY  OTHER  CAR  ) 


3-42 


09/11  /74 


AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MOOEL 


CHART  TITLE  - SUBROUTINE  INPUT 


/ INPUT  / 


02.05- 


->* 


R E AOS  TRAIN  COLLISION 
MOOEL  INPUT  FROM 
CAROS  ANO  PRINTS 
INPUT 


READ  INPUT 


RUN  TITLE  CARO.  80 
CHARACTERS 


I 


01 


/ READ  FROM  DEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 99  / 

/ INTO  THE  LIST  / 


I 


LIST  = 

( TITLE( I ) , I = 

1 .20) 


PROGRAM  CONTROL 

I 

I 


/ RE AO  FROM  DEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 98  / 

/ INTO  THE  LIST  / 


I NOTE  04 

• LIST  = TO,  TF,  ♦ 

• DT,  DTP,  MODE  ♦ 

• •***•*•*•• 

I 

I 

I 05 

/ READ  FROM  DEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 97  / 

/ INTO  THE  LIST  / 

I 

I 

| NOTE  06 

*****••**•• 

* LIST  = NCARS,  • 

* WEIGHT,  VW,  VTRN,  * 

♦ ACO,  APO  • 

♦ *•**•**•*• 

I 


I 


07 


/ READ  FROM  DEV  / 
/ 5 / 

/ VIA  FORMAT  / 
f 96  / 

INTO  THE  LIST  / 

I 

I 

I NOTE  08 

LIST  = • 

CSPACd  ),I  = 1,6)  * 


' 1 

09  | 

1 NOTE  18 

/ READ  FROM 

DEV  / 1 

• LIST  = * 

/ 5 

/ 1 

• (DEL  ( 1 , 2 ) ,1  = * 

/ VIA  FORMAT 

/ | 

* 1,10)  * 

/ 96 

/ 1 

/ INTO  THE  LIST 

/ | 

1 

1 

1 

1 19 

1 NOTE  10  | 

/ READ  FR CM  DEV  / 

******** 

* * • | 

/ 5 / 

* LIST  = (CRU)  , 

I - • | 

/ VIA  FORMAT  / 

• 1,5) 

• | 

/ 95  / 

i 

1 

/ INTO  THE  LIST  / 

EACH  FCO  AND  DEL  | 

1 

I NOTE  20 

REQUIRES  2 CARDS 

WITH  | 

2 VALUES  ON  SECOND  1 

* LIST  = UR  AT  E(  1),  * 

CARD 

1 

• UR AT  E ( 2 ) * 

BLANK  CARD  MUST 

BE  I 

FURNI SHED  IF  NO 

DATA  | 

1 

FOR  LAST  TWO  ENTRIES.  1 

1 

i 

i 

IF  INDICATES  WHICH 

1 

11  1 

FCO-OEL  TABLE  USED 

— 

| 

FOR  EACH  SPRING. 

/ READ  FROM 

DEV  / | 

MUST 

/ 5 

/ | 

BE  1 OR  2 

/ VIA  FORMAT 

/ 1 

1 

/ 95 

/ 1 

| 

/ INTO  THE  LIST 

/ | 

1 21 

| 

1 

/ READ  FR  CM  OEV  / 

| 

1 

/ 5 / 

1 NOTE  12  1 

/ VIA  FORMAT  / 

/ 94  / 

• LIST  = 

♦ 1 

/ INTO  THE  LIST  / 

• ( FC  0 ( I , 1 ) , I 

— 

* 1,10) 

• i 

1 

1 

1 

I NOTE  22 

13  1 

* LI  ST  = ( I F ( I ) , I = * 

— 

| 

* 1, NCARS)  * 

/ READ  FROM 

DEV  / | 

/ 5 

/ 1 

1 

/ VIA  FORMAT 

/ 1 

1 

/ 95 

/ 1 

| 

/ INTO  THE  LIST 

/ 1 

WFACTR  READ  ONLY  IF 

1 

1 

1 

MODE  = 1 

1 

1 1 
I NOTE  14  1 

] 

* * * • * ♦ * * 

♦ * • | 

• 23 

• LIST  = 

• 1 

* * 

* ( FC  0 ( I , 2 ) , I 

* • 

* 1,10) 

• i 

FALSE  * * 

♦- 

* MODE  .EQ.  1 » 

1 

l 

1 

* * 

1 

i 

1 

* * 

1 

15  | 

1 

* * 

— 

i 

1 

• 

/ READ  FROM 

DEV  / 1 

1 

1 TRUE 

/ 5 

/ 1 

1 

1 

/ VIA  FORMAT 

/ 1 

1 

i 

/ 95 

/ 1 

1 

1 

/ INTO  THE  LIST 

/ 1 

1 

1 

i 

1 

1 

1 24 

I NOTE  16  1 

1 

/ READ  FR  CM  DEV  / 

******** 

* * * | 

1 

/ 5 / 

* LIST  = 

* 1 

1 

/ VIA  FORMAT  / 

• (DEL (1,1), I 

1 

/ 95  / 

* 1,10) 

* 1 

1 

/ INTO  THE  LIST  / 

«*•***«* 

1 

* * * | 
1 

1 

1 

1 

1 

1 

1 

17  | 

1 

1 

1 

I NOTE  25 

| 

1 

*•*****♦*•» 

/ READ  FROM 

DEV  / 1 

i 

• LIST  = WFACTR  * 

/ 5 

/ 1 

i 

/ VIA  FORMAT 

/ 1 

1 

1 

/ 95 

/ 1 

♦- 

> | 

/ INTO  THE  LIST 

/ 1 

1 

— 

1 

PRINT  INPUT 
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AUTOFLOW  CHART  SET  C-wc* 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - SUBROUTINE  INPUT 


13.26 >♦ 

I 01 

/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1000  / 

/ FROM  THE  LIST  / 


I 

I 

I NOTE  02 

*********** 

* LIST  = * 

* (TITLE!  I ) »I  «=  • 

* 1,20)  * 

*********** 

I 

I 03 

/ WRITE  TO  DEV  / 

/ 6 / 

/ VIA  FORMAT  / 

/ 1001  / 

/ FROM  THE  LIST  / 


I 

I 

I NOTE  04 


* LIST  = TO,  NCARS,  * 

* TF,  WEIGHT,  DT,  * 

* VW,  DTP,  VTRN,  * 

* MODE,  ACO,  APO  * 

*********** 

I 

I 

I 05 


/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1002  / 

/ FROM  THE  LIST  / 


I 

| NOTE  06 

*********** 

* LIST  = * 

* ( FCO ( I , 1 ) , DEL ( I,  * 

* 1),FC0(I,2),  * 

* DE  L ( 1 , 2 ) « 1 = * 

* 1,10)  * 

*********** 

I 

I 

I 07 


/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1003  / 

/ FROM  THE  LIST  / 


I 

I 

| NOTE  08 

* ********** 

* LIST  = URATEd  ),  * 

* URATE ( 2 I * 

*********** 

I 

I 09 

/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1004  / 

/ FROM  THE  LIST  / 


< >* 

I 111 


/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1005  / 

/ FROM  THE  LIST  / 


I 

I 

I NOTE  12 

*•*****•••* 

* LIST  = ( IF ( I ) , I = * 

* 1 , NCARS I * 

*********** 

I 

I 

I 13 


/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1006  / 

/ FROM  THE  LIST  / 


I 

I 

I NOTE  14 
*****••*••* 

* LIST  «=  * 

* ( SPAC  ( 1 ) ,CR  ( I ) , * 

* I = 1,5),  SPAC (6 ) * 
*********** 

I 

I 

* 15 

* * 

• * 

FALSE  * * 

* MODE  .EQ.  1 * 

* * 

* * 

* * 

* 

ITRUE 


I 

I 

I NOTE  10 

*********** 

* LIST  * (i,i  « * 

* l, NCARS)  * 

*********** 

I 


16 


/ WRITE  TO  DEV  / 
/ 6 / 

/ VIA  FORMAT  / 

/ 1007  / 

/ FROM  THE  LIST  / 


I 

I 

I NOTE  17 

*********** 
* LIST  = WFACTR  • 
*********** 
I 

I 18 
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CHART  TITLE  - NON-PROCEDURAL  STATEMENTS 


99 

90 

97 

96 

93 

9*» 

1000 

1007 

1001 


1002 


1003 

1004 
1003 
1006 


COMMON  IMP (20) , WOT (20 ) ,FC (2 1 ) , XC ( 20 ) , VC ( 20 ) * AC ( 20) * X P ( 20 » , VP ( 20 ) , 
AP120),VIMP(20,6),DVP(20)»DP(20)*DM(20)*FM(20)*URATE(2)» 

OR (20) ,DELI10,2),FC0( 10*2) ,IF(20) .TITLE (20) ,TO,NC ARS.TF, 
WEIGHT,DT,VW,DTP,VTRN,M00E,AC0,AP0,SPAC(6 ) ,CR (5) • WFACTR, 

XW  ,T,G,N,AIMP(20,6,5> ,SI (20,6,5) 

FORM  A T ( 20A4  ) 

FORM  AT(4F10.0« 15) 

FORMAT  I I5,5X,5F10.0) 

FORM  ATI 6F10  *0) 

FORMAT! 8G10.0) 

FORMAT! 2014 ) 

FORMAT! lHlt  40X  * 2 1HTR  A IN  COLLISION  MODEL  / 20X.20A4) 

FORM  ATI 1H0,  //  9H  WFACTR  =,F10.4) 

FORMAT! 1H0, 25X  » 4HT0  = ,F8.4,4H  S EC , 8 X , 1 3HN0.  OF  CARS  = ,I4,  / 

26X  , **HTF  =tF8.4,4H  SEC  ,9  X , 1 1 HCAR  WE  I GHT  *,FI2.1*4H  LBS  / 
26X.4HDT  «,F8.5,4H  SE C .3 X , 1 8HBARR I ER  VELOCITY  =*,Fi2.2,4H  MPH/ 
20X.10HDT  PRINT  =,F0.4,4H  S EC » 5X  , 16HTR A IN  VELOCITY  =,F12.2, 

4H  MPH  / 

24X,6HM00E  = , 14, 15X , 14H TR A IN  DECEL.  «,F12.2»4H  G#S  / 

45X, 18HPASSENGER  DECEL.  =,F12.3,4H  G»S  ) 

FORMAT! 1H0,  //  10X,38HTRA1N  FORCE-DEFLECTION  CHARACTERISTICS  / 
I3X,6HFC0( 1 ) « lOXt  6HDEL 1 1 I , I OX ,6HFC0 ( 2 ) * 10X , 6HDE L( 2 I / 

14X,3HLBS,13X,3HFT.,13X,3HLBS,13X,3HFT.  // 

I 10X,F10.1,8X,F7.2 ,7X,F10.1 ,8X,F7.2>  ) 

FORMAT! 1H0, 10X.26HUNL0ADING  RATES:  URATE (1) « * F 14 .2 ,5HL 8/ FT  , 
1GX,9HURATE(2)«,F14.2,5HLB/FT  /) 

FORM  AT ( 1H0*  //  5X.7HCAR  N0.,28X,15I5  ) 

FORMAT ( 1H  ,4X,35HF0RC  E-OEFLECT ION  CHARACTERISTIC  NO.  ,1515  ) 

FORMAT! IHO,  //  IOX , 1 7H PASSENGER  SPAC ING , 13X , l 6HP AS SE NGER  DECEL.  / 
1 OX, 1 3HD 1ST  ANC  E - FT, I7X , 13HDI ST ANCE  - IN  / 

( 13X,F0.2,22X,F8.2  ) ) 
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AUTOFLOW  CHART  SET 


TRAIN  COLLISION  MODEL 


CHART  TITLE  - SUBROUTINE  CLEAR ( A ,B  ) 


/ CLEAR  / 


02.01 >* 


CLEARS  (SETS  TO  ZERO! 
A BLOCK  OF  STORAGE 
IDENTIFIED  BY 
THE  ADDRESSES  OF  THE 
TWO  ARGUMENTS 
CALL  CLEAR ( P * Q ) WILL 
CAUSE  ALL  BYTES  TO  BE 
SET  TO 

ZERO  FROM  ADDRESS  P 
THROUGH  THE  FULL-WORD 
ADDRESS  AT  0 


* * FALSE 

5(1)  .EQ.  0.0  * 


I I = I ♦ 1 | 

I I 

| A ( I ) = 0.0  | 
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CHART  TITLE  - NON-PROCEDURAL  STATEMENTS 


DIMENSION  AllliBIll 


3-47/3-48 


DOT  LIBRARY 


III 

III 

II 

II 

III 

D 

]Q[ 

7/ 

